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Abstract 

We discuss finding large integer solutions of d = 2x 3 + y 3 + z 3 by using 
Elsenhans and Jahnel's adaptation of Elkies' LLL-reduction method. 
We find 28 first solutions for \d\ < 10000. 



1 Introduction 

Fermat's Last Theorem states that 

x 3 + y 3 = z 3 
has no solutions in non-zero integers. 
Over the years, there has been interest in small values of 

d = x 3 + y 3 - z 3 = x 3 + y 3 + (-z) 3 (1.1) 



Since u 3 = 0, 1,8 (mod 9) for all integers u, it is an easy enumeration of 
possibilities to see that NO solutions can exist for d = 4, 5(mod9). It is 
conjectured that solutions exist for all other values of d and it is an on-going 
project to find such solutions for small d. 
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For 



d = 2x 3 + y 3 + z 3 



(1.2) 



there are no congruence restrictions on d. The problem is discussed in Guy 
[I], and Koyama presents a specialized search method in [5] but much less 
work seems to have been done on this problem. It appears somewhat easier 
than (1.1) as there are only 3 values of \d\ < 1000 for which no solution is 
known, namely d = 148, 671, 788, and only a total of 62 values for \d\ < 10000, 
according to Hisanori Mishima's web-site [6]. 

Recent, Elsenhans and Jahnel [3] used Elkies' method [2] to discover previ- 
ously unknown first solutions to (1.1) for some \d\ < 1000, but with x, y and 
z very large in size. The paper [3] does not contain many technical details. 
Jahnel's web-site gives a zipped file containing the C-code used, together 
with a small Word document describing the method. Unfortunately, all the 
comments in the code and most of the Word file are in German, which the 
author last seriously studied in the 1970s. 

In an attempt to reverse-engineer their methodology, equation (1.2) was used 
as a test-case. Since this problem might be popular with amateur mathe- 
maticians, the following description is intended to be as simple as possible. 



2 Basic Theory 

We are considering solutions which cannot be found by simple searching over 
small values of x,y,z. We assume \d\ < 10000, and min(|y|, \z\) > 100, so 
that there must be a mixture of positive and negative values of x, y, z. 

We write the problem as 

2x 3 + y 3 = z 3 + d (2.1) 

with d positive or negative, and, thus, we can assume, without loss of gener- 
ality, that z > 0. Then at least one of x, y must be positive. 



:) +2(-) =1 + ^ (2-2) 



We have, 

:!)-© 3 -+p 

and, if z is large, we are effectively looking for rational points close to the 
curve 

Y 3 = 1 - 2X 3 (2.3) 
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and it is for this type of problem that Elkies initially proposed his method. 
We have 3 intervals of interest 

(a) If X > 0, Y > 0, then < X < ^0T5. 

(b) If Y < 0, then X > ^53. 

(c) For any X < we have Y > 0. 

Define 

Y = f(X) = (1 - 2X 3 ) 1/3 (2.4) 
and consider, for ease of exposition, < X < v^OTo" « 0.7937. 

We have, 

dY _ -2X 2 d 2 Y _ -AX 

dX ~ (1 - 2X 3 ) 2 / 3 ~ (1 - 2X 3 ) 5 /3 

Let X be a random point in the interval (0,0.7937), with Y = f(X ) so 
that 

dY -2X1 d2Y - 4X o 

dX ~ F 2 ~ F 5 

Consider the interval / = [X — h/2, X + h/2] with h strictly positive, where 
we usually assume the end-points stay inside the basic interval for X. In this 
interval, we wish to surround the curve by a parallelogram. Elsenhans and 
Jahnel use the German word fliesen = flagstone for such a parallelogram. We 
assume the parallelogram has 2 edges parallel to the F-axis, but the other 2 
edges are angled downwards. In fact, we only need to specify the equation 
of the middle downward line of the parallelogram. 

We know the gradient of the curve at (X , Y ) is A = — 2Xq/F 2 , so we look 
for sides of the parallelogram with gradient A. The equation of the tangent 
itself is 

Y = AX + (Y -AX ) (2.5) 
but the properties of the curve show that this line is always above the curve. 

What we want is a line 

Y = AX + B 

with B < Y — A X which lies both above and below the curve in /. 
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We have, using Taylor series, 



which is just 



f(X + \) « f(X ) + + /"(*o)y 



A (X + h/2) + (Y -AX + h 2 f"(X )/8) 



(2.6) 



If we set B = Y - A X + h 2 f"(X ) /8, the resulting line lies totally below the 
curve, remembering that the second derivative is negative, but if we define 

B = Y -AX + h 2 f"(X )/lQ 

we get a line which satisfies the requirement of above and below the curve. 

Using the above discussion, we look for rational X = x/z,Y = y / z with 

\X-X \<h/2 (2.7) 

and 



\Y-(AX + B)\ < K 
where K > but usually very small. 

If we restrict z to < z < L, this gives 

\x-X z\ < h\z\/2 < hL/2 

and 

\y - Ax - Bz\ < K\z\ < KL 
We, thus, have the three restrictions 



2 2X 
hL X ~ ~hL Z 



< 1 



and 



and 



A 

-y - ~77t x - 



B 



KL KL KL 



< - < 1 

Lj 



< 1 



(2.8) 



(2.9) 
(2.10) 



Define the matrix 



F = 



( 


2 

hL 





-2X 
hL 


\ 




-A 


1 


-B 






KL 


KL 


KL 




V 








1 
L 


/ 



with the columns denoted / l , / and / . 
Then we are looking for vectors in the lattice 

v = xf x + yf_ 2 + zf_ z x,y,z G Z 
which have small components. 

This is what the LLL algorithm does, see Cohen [TJ. So we apply LLL- 
reduction to F, and, hopefully, find a solution to the problem with \d\ small. 

As an example, let h = 0.001, K = 0.00001 and L = 1000. X = 0.31415 
gives A = -0.205986 and B = 1.043599, so that we have 

/ 2 -0.6283 \ 

F= 20.5986 100 -104.3599 
\ 0.001 / 



Applying the function qflll from the software package Pari version 2.3.4 gives 
the integer matrix 

/ -15 -74 -313 

M= -47 -230 -976 

\ -48 -235 -997 

and we have 

0.1584 -0.3495 0.4151 

/•' U = | 0.2962 0.2801 -0.5415 

-0.0480 -0.2350 -0.9970 

Investigating the first column of M quickly gives 

2(-15) 3 + (-47) 3 + 48 3 = 19 
whilst the second column gives 

2(-74) 3 + (-230) 3 + 235 3 = 427 

For the other two intervals X > v^O.5 and X < 0, the only difference is the 
behavior of the derivatives, but the resulting formulae all turn out the same. 
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3 Pyramid Searching 



The first version of the algorithm, which we implemented in Pari, generated 
F, computed M and just investigated the 3 columns separately as possible 
solutions. The columns, however, give small vectors in the lattice but not 
necessarily linked to small values of \d\. 

To ensure we do not miss such solutions we need to investigate further. As 
pointed out by Elsenhans and Jahnel, the relations (2.9) and (2.10) give a 
feasible region which is, in fact, a pyramid. The vertex is at the origin, and 
the four other vertices at the solutions of 















H 













(3.1) 



A look at the numerical solution in the previous sections shows that there can 
be a huge number of lattice points within this pyramid, so that testing them 
all will be impractical. The matrix M, however, is unimodular (determinant 
= ±1), so 




with x', y', z' G Z. The matrix H = F M is much nicer to deal with and gives 
a much smaller pyramid to search. 

We work out the equations of the 8 edges of the pyramid, and loop through 
the possible z' values. For each z' =constant plane, we work out the (x', y') 
limits and loop through these. This does not work for z' = so we specify 
the (x',y', 0) triples. We, especially, do not want to miss the first column of 
M as it is, very often, the smallest non-zero vector in the lattice. 



4 Practical Implementation 

If we study the numerical results in Koyama, we find that the solutions all 
have X roughly in the range [—10, 10], so we generate random values of X in 
such a range. For input values of Xq, h,K,L a Pari 2.3.4. program generates 
F and the Pari code qflll generates the LLL-reduction matrix M and, hence, 
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F M. We then use the pyramid method of section 3 to search for a solution 
amongst the 62 unsolved values of d from Mishima's web pages. 

A large amount of experimental testing suggests that K = 0(h 2 ) and L = 
0(1/ h) are good choices, but these are fairly imprecise relationships. 

The program is very small and was run on several lap-tops. The results found 
so far are 



TABLE 1 
Solutions of d — 2x 3 + y 3 + z 3 
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